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A phenomenological model defined in a two dimensional lattice is employed to investigate the 
d-wave superconductivity and the competing antiferromagnetic order in cuprate superconductors. 
Near the optimally doped regime, we show that it is possible to induce the spin density wave (SDW) 
and the accompanying charge density wave (CDW) orders with stripe modulations by applying a 
magnetic field. The periods of the magnetic field induced SDW and CDW are 8a and 4a, respectively. 
The spatial profiles of the local density of states (LDOS) at various bias energies have also been 
numerically studied. Near and beyond the energies of the vortex core states, we found that the 
LDOS may display the CDW stripe-like modulation throughout the whole magnetic unit cell. For 
energies closer to the zero bias, the stripes appear to be rather localized to the vortex. The intensity 
of the integrated spectrum of the LDOS shows that the strength of the stripe modulation is decaying 
away from the vortex core. This feature is in good agreement with STM experiments. The case for 
the magnetic-field induced SDW/CDW with 4- fold symmetry has also been studied. 
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In hole-doped cuprate superconductors, the interplay 
between the d-wave superconductivity (dSC) and anti- 
ferromagnetic (AF) order has been studied extensively in 
the literatures. Inelastic neutron scattering experiments 
showed the presence of incommensurate spin structures 
in Lao- zSr^CuCU (LSCO) Q with a spatial periodicity 
8a in the presence of a magnetic field. In addition, NMR 
imaging experiment on YBa2Cu307_ :E (YBCO) observed 
a strong AF fluctuation outside the vortex Q indicating 
the possible existence of spin-density-wave (SDW) gap 
inside the vortex. A recent STM experiment by Hoff- 
man et al. |j| seems to confirm the coexistence of the 
static charge modulation and the superconductivity on 
BiaSraCaCuaOg (BSCCO) under a magnetic field. The 
authors reported that a four unit cell checkerboard pat- 
tern is localized in a small region around the vortex, and 
its intensity is exponentially decaying away from the vor- 
tex core. All these experiments indicate that the AF 
fluctuations could be pinned by the vortex cores and they 
may form a static SDW/chargc-dcnsity-wave (CDW) like 
modulations in certain samples of cuprate superconduc- 
tors. Theoretically, a number of works proposed that the 
observed checkerboard patterns could be explained by 
the SDW order with the two-dimensional (2D) 0, H Q 



or stripe pUSlISi modulations. Both the 2D- and stripe- 
SDW orders induced by the magnetic field are well known 
to have the accompanying CDW modulations. In par- 
ticularly, the checkerboard pattern has been attributed 
to the superposition of stripe modulations of the CDW 
0, E3 oriented along x- and y- directions. In all these 
studies, the STM spectra obtained from the experiment 
have been directly interpreted in terms of the CDW order 
induced by the magnetic field. This is because that the 
conventional wisdom leads us to believe that the symme- 
try of the CDW order should be reflected in the STM 



spectra. On the other hand, the CDW order represents 
the charge density configuration in the ground state while 
the STM spectrum or the spatial distribution of the local 
density of states (LDOS) is determined by the behavior of 
the low-energy excitations in the system. It is therefore 
extremely interesting to know the difference between the 
CDW order and the spatial profile of the LDOS. This 
issue has not yet been addressed in the literature, and 
thus is the main purpose of the present paper. 

The phonemenological t — t' — U — V model defined 
in a 2D lattice and the Bogoliubov-de Gennes' equations 
will be used to numerically study the interplay between 
the d-wave superconductivity (dSC) and the competing 
AF order for samples close to the optimal doping. First 
the phase diagrams for the coexistence between these two 
orders as a function of U/V in both zero and finite mag- 
netic field are going to be examined. Two different values 
of U /V will be chosen such that their magnitudes yield 
only the dSC when the magnetic field is zero. In the pres- 
ence of a magnetic field, the U/V(= 2.39) will give rise 
to a 2D SDW/CDW-like modulations while a stronger 
U/V(= 2.44) would make the induced SDW/CDW to 
have one dimensional stripe- like structures. Then the 
spatial distributions of the LDOS as a function at vari- 
ous bias energies will be calculated and compared with 
those of the CDW. The spatial distributions of the inte- 
grated LDOS have also been obtained, and the result for 
U /V = 2.39 shows strong LDOS intensity near the vor- 
tex core and its distribution seems to be round, not the 
four-fold symmetry as expected for the CDW and dSC 
orders. The result for U/V = 2.44 shows that the stripe- 
like modulations are still existing but they are more lo- 
calized near the vortex core and this is different from 
the CDW where the stripes are extended over the whole 
magnetic unit cell. These features are in good agreement 
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with the experiments of Pan et al. ^lj an d Hoffman et al. 

a 

We start with an effective mean- field t — t' — U—V 
Hamiltonian in the mixed state by assuming that the 
on-site repulsion U is responsible for the competing an- 
tiferromagnetism and the nearest-neighbor attraction V 
causes the d-wave superconducting pairing. 

H = - tij c| CT c jCT +^2(U (n is ) - M)4 CT c iCT 

ijer ier 

+ ^(A ijC t T c] i + Ar j c ji c iT ) , (1) 



where iy is the hopping integral, /i is the chemical po- 

CijCj|) is the spin-singlet 



tential, and A; 



Tl c iT c ji 



d-wave bond order parameter. The Hamiltonian above 
shall be diagonalized by using Bogoliubov-de Gennes' 
(BdG) equations, 
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where Wye, = -iy + (U(n ia ) -^)<5y. Here, iy = (iy)e t¥3i j . 
The Peierl's phase factor t^y = ^ J r rj A(r) • dr , with 
the superconducting flux quantum <J>o = hc/2e. Within 
the Landau gauge A(r) = (—By, 0,0), each magnetic 
unit cell can accommodate two superconducting vortices. 
The vortex carries a flux quantum <I>o and locate at the 
center of a square area containing N x /2 x N y sites with 
N x = 2x N y . Here, we choose the nearest-neighbor hop- 
ping (iy) = t = 1 and the next-nearest-neighbor hop- 
ping (iy) = t' = —0.25 to match the curvature of the 
Fermi surface for most cuprate superconductors. The 
exact diagonalization method to self-consistcntly solve 
BdG equations with the periodic boundary conditions 
is employed to get the N positive eigenvalues (E n ) with 
eigenvectors (it£j-, vQ) and negative eigenvalues (E n ) with 
eigenvectors (— u-l*, u™?). The self-consistent conditions 
are 
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Ay = T W + V ?X) tanh(^) , (3) 

n=l 

where u" = (— uV\) and v™ = (u™£,v^) are the row 
vectors, and f(E) = l/(e l3E + 1) is Fermi-Dirac dis- 
tribution function. Since the calculation is performed 
near the optimally doped regime, the filling factor, rif = 
Si(r( c L c io-} /N X N V , is fixed to be 0.85, i.e., the hole dop- 
ing 5 = 0.15. Each time when the on-site repulsion U is 
varied, the chemical potential [i needs to be adjusted . 

In the limit of U/V <C 1, the system is in the state 
of pure dSC. In the opposite limit U/V 3> 1, the sys- 
tem is the SDW states. However, the most anomalous 
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FIG. 1: The maximum value of the staggered magnetization 
in the zero-field (a), and in the finite field (b). In (b), the open 
circle and the solid circle represent the (maximum) staggered 
magnetization at the center of the vortex core and away from 
the vortex cores, respectively. Stripe^) is the region where 
the stripe modulations are intrinsic and stripej//) corresponds 
to the region where the stripe modulations are field-induced. 
The size of the unit cell is N x x N y — 48 x 24 correspoding 
to a magnetic field B ~ 24T. 



properties of cuprate superconductors do not correspond 
to these extreme limits, but are in the intermediate case 
where both the SDW and the SC may coexist. In or- 
der to simplify our discussion of the ratio of U to V, we 
set V = 1.0. In Fig. 1(a), without the magnetic field, 
B = 0, the staggered magnetization (M\) or the SDW 
shows either the stripe modulation or the uniform distri- 
bution with Mi = in the background of dSC depending 
upon the magnitude of U. For U > U c \, the staggered 
magnetization exhibits the stripe modulation with 8a as 
its periodicity (a is the lattice constant). On the other 
hand, for the U less than U c i, it shows the uniform dis- 
tribution which is equivalent to the state of pure dSC. It 
is important to point here that a two-dimensional (2D) 
SDW modulation can never be obtained when B = in 
the present self-consistent calculation. The transition be- 
tween the SDW-stripe modulation and the uniform dis- 
tribution is discontinuous. In Fig. 1(b), under an ap- 
plied magnetic field, the staggered magnetization (solid 
line) displays the stripe modulation, the two dimensional 
SDW, or the uniform distribution depending on U. As U 
is greater than U c \, the stripe modulation existed in the 
zero field would be slightly enhanced under a magnetic 
field (Stripem). In the region of U C 2 < U < U c \, the field 
induced staggered magnetization shows a stripe modula- 
tion which disappears in the zero field (Stripe^/)). For 
B / 0, the transition between stripem and stripe(m is 
not apparent and only the slope shows a weak discontinu- 
ity. When U < U C 2, the field induced AF order changes 
from the stripe-like to a 2D SDW. If U goes down far 
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below 2.3, the AF order could be completely suppressed 
both inside or outside the vortex cores. The transition 
between the stripe modulation and the 2D SDW is of the 
first order. Accompanying the stripe-like AF order, there 
also exists a CDW with the stripe modulation of the pe- 
riod 4a. At the same time the dSC order parameter also 
acquires a CDW-like stripe modulation. In Fig. 1(b), 
the staggered magnetization (open circles) at the vortex 
core center seems to be weakly U dependent. In the fol- 
lowing, U = U a = 2.44 and U = Ub = 2.39 are chosen 
for our study of the LDOS. Both of these U values would 
not generate the SDW/CDW order for nearly optimally 
doped cuprate superconductors when B = 0. In order to 
understand the characteristics between the two cases, we 
start with the LDOS formula 

1 2N 2 



n .k 



n , k 



/'(E n , k + E)} , 



(4) 



where pi(E) is proportional to the local differential tun- 
neling conductance as measured by STM experiment, and 
the summation is averaged over M x x M y wavevectors in 
first Brillouin Zone. 

The LDOS as a function of energy have been numeri- 
cally calculated at the vortex core center (solid line) and 
at site far away from the vortex (dashed line) for Ua and 
Ub, the results are respectively given in Fig. 2(a) and 
2(b). There the spatial profiles of the field induced CDW 
modulations arc presented in Fig. 2. For Ua, the field 
induced CDW has the stripe-like structure which extends 
over the whole magnetic unit cell with a period 4a while 
for Ub, the CDW modulation becomes two dimensional 
with four-fold symmetry. 

Since the SDW gap develops inside the vortex core, the 
resonance peak appeared in the LDOS ^3 near the zero 
bias in a pure dSC at the core center is suppressed and 
splits into two peaks [l4| . This feature can be seen in the 
Fig. 2. Furthermore, STM experiments on YBCO [13 
and BSCCO [ill measured the double peaks structure 
within the maximum superconducting gap. The states 
associated with these two peaks have been referred as 
the vortex core states. Here we would like to point out 
that our results for the LDOS seem to agree better with 
YBCO [13 than BSCCO In Fig. 2, those vortex 

core states with negative energies are centered at Ea = 
-0.19 with a width AE A = 0.08 and E B = -0.21 with 
AEb = 0.04 for Ua and Ub, respectively. 

Because the energy dependence of the LDOS at the 
vortex center does not make any clear distinction be- 
tween Ua and Ub, we examine the spatial profile of the 
LDOS at various energies and look for their differences. 
In Fig. 3, the LDOS maps (the spatial distribution of the 
LDOS), which arc the LDOS with fixed energy at each 
site of the magnetic unit cell, have been calculated at 
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FIG. 2: The LDOS as a function of energy, left panel, at the 
vortex core center for (a) Ua = 2.44, (b) Ub = 2.39. The 
solid line is at the vortex core center, and the dash line is at 
the site far away the vortex. The arrow points the vortex core 
states (VCS). The spatial profiles of the corresponding CDW 
are shown in the right panels of Fig. 2. The wavevectors in 
first Brillouin Zone are M x x M v = 24 x 24. 



energies ranging from 0.0 to —0.4 with SE = 0.01 decre- 
ment. For Ua, as the energy far below the vortex core 
states (VCS) and close to the zero bias, the pattern in 
Fig. 3(a) shows that the stripe structure with periodic- 
ity 4a is strongly localized near the vortex core, and its 
strength drops dramatically away from the vortex. When 
the energy near the VCS, such as Fig. 3(b), the stripe 
modulation in the LDOS extends over the whole unit 
cell, similar to the feature in CDW in Fig. 2(a). If the 
bias voltage goes above or becomes more negative than 
the VCS, such as Fig. 3(c), the stripe modulation is still 
extensive, but the intensities inside the vortex core are 
depressed than those outside the vortex core. Thus the 
STM spectra or the LDOS do not have all the features 
of the CDW. 

For Ub , as the energy far below the vortex core states 
(VCS) [see Fig. 3(e)], the LDOS pattern shows a round 
bump with the size of a vortex core. When the energy is 
near the VCS, such as in Fig. 3(f), the LDOS shows a dis- 
tribution with rather weak oscillations, and its feature is 
not quite similar to that of the CDW in Fig. 2(b). If the 
energy is above the VCS, such as that in Fig. 3(g), the 
modulation here clearly displays the pattern with 4-fold 
symmetry. However, the STM images obtained experi- 
mentally [3 are results of integrating the spectral density 
between energies E± and Ei, which is defined as 



S(E U E 2 ) 



E 2 

J2pi(E)6E, 



(5) 



There, |3j E\ is taken to be and E<i has been set near 
the energy of the VCS below the chemical potential. In 
Fig. 3(d), the integrated spectrum 5(0.0, —0.23) of the 
LDOS is obtained by summing the LODS from E\ to 
Ei with an energy spacing 0.01. It shows that the inten- 
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FIG. 3: The LDOS maps are at energies (a) and (e)-close 
to the zero bias, (b) and (f)-near the VCS, and (c) and (g) 
above the VCS. The integrated spectrum of LDOS S(Ei,E 2 ) 
of (d) and (h) are from the Fermi energy to the upper bound 
of the VCS, i.e., 5(0.0, -0.23). The left and right panel are 
for Ua = 2.44 and Ub = 2.39, respectively. The wavevectors 
in first Brillouin Zone are M x x M y =6x6. 



as expected for a dSC [TJ|. This feature is also consistent 
with the STM experimental measurements provided 
that the samples used there is different from that of f|. 

In conclusion, we have numerically investigated the in- 
terplay between the dSC and the competing SDW/CDW 
orders by varying the strength of U /V for samples close 
to the optimal doping. For finite field, we show that both 
stripe-like and two dimensional SDW/CDW orders, de- 
pending on the magnitude of U /V, may be induced in the 
background of the dSC. We in particular calculate the 
spatial distribution of the LDOS with and without the 
SDW/CDW stripes at various bias energies. We point 
out that the checkerboard pattern near the vortex core 
observed by Hoffman et al. 0] could be interpreted in 
terms of superposition of the field-induced x- and y- ori- 
ented stripes. The obtained features of our integrated 
LDOS with stripes and without stripes [see Figs. 3(d) 
and 3(hV| are in good agreement with the STM exper- 
iments yl [n|. It is well known that the theoretically 
obtained profiles for the dSC order parameters together 
with the induced 2D SDW/CDW modulations near the 
vortex core exhibit a clear 4-fold symmetry. This sym- 
metry, however, so far has not been confirmed by ex- 
isting STM measurements. Here, we predict that this 
4-fold symmetry should be more easily detectable when 
the bias is placed beyond the energies of the vortex core 
states [see Fig. 3(f)]. 
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sity of the stripe modulations is dominantly concentrated 
near the vortex core and decays rapidly when away from 
the vortex. This feature originates from the behaviors 
of the LDOS at energies with magnitudes smaller than 
those of the VCS [see Fig. 3(a)]. Since the LDOS maps 
display the properties of the eigen function at the energy 
E, the checkerboard pattern observed by the experiments 
could be explained in terms of the superposition of the 
degenerate eigenfunctions describing the stripe modula- 
tions along x- and y- directions. All these behaviors are 
in good agreement with the experiment Q . We expected 
that the qualitative features obtained for Ua should still 
remain even if U > U c i and the stripe phase is intrinsic 
not magnetic field induced (see Fig. 1). In Fig. 3(h), 
the integrated spectrum of the LDOS for U b exhibits a 
round bump over the vortex core region which is differ- 
ent from that of the CDW as shown in Fig. 2(b). For a 
sample without the stripe modulations, the profile of the 
integrated LDOS near the vortex core seems to be rather 
round and does not possess the strong 4-fold symmetry 
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